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Abstract 

In this paper we provide a closed mathematical formulation of our previous results 
in the field of symbolic dynamics of unimodal maps. This being the case, we dis- 
cuss the classical theory of applied symbolic dynamics for unimodal maps and its 
reinterpretation using Gray codes. This connection was previously emphasized but 
no explicit mathematical proof was provided. The work described in this paper not 
only contributes to the integration of the different interpretations of symbolic dy- 
namics of unimodal maps, it also points out some inaccuracies that exist in previous 
works. 
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1 Introduction 



A symbolic sequence is a transformation of a sequence of real numbers into a 
sequence consisting of a set of symbols. Regarding unimodal maps, the car- 
dinality of that set is two and it is determined by the turning point of the 
iteration function of the map. Accordingly, each symbol represents the relative 
position of a real- value with respect to the turning point. In [11] it is pointed 
out the existence of an inner order of the symbolic sequences, along with the 
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relationship between this order and the initial condition and the control pa- 
rameter of the underlying chaotic system. The considerations and results of 
[11] were later improved and enlarged through different contributions, being 
the most important [8] and [13]. In [ i] it was remarked that the order of the 
symbolic sequences can be interpreted using the concept of Gray codes. In this 
novel approach to the problem, the symbolic sequences are finally converted 
into a figure which is a real number between and 1 called Gray Ordering 
Number or simply GON. Afterwards , [9] drew the bridge between the ideas 
of [4] and the main theory of applied symbolic dynamics as expressed in [13]. 
Finally, some theorems are offered in [14], which enlarge the theoretical frame- 
work of the GON of unimodal maps. In [ ] it is explained that the dynamical 
properties of unimodal maps by means of the GON are a translation of the 
theoretical framework inherited from [11]. Nevertheless, there is no direct and 
explicit proof of this equivalence. One of the main applications of the concept 
of the GON is the estimation of the control parameter of unimodal maps for 
cryptanalysis [2,7,12,6]. The precise definition of the key space of a cryptosys- 
tem is a commitment in cryptography. In the context of chaotic cryptography, 
it implies that the control parameters and initial conditions of the chaotic sys- 
tem must be selected to guarantee chaoticity, and to avoid the estimation of 
either control parameters or initial conditions from partial information about 
the chaotic orbits [1, Rule 5]. In case that this partial information arises from 
the symbolic sequences of the chaotic map used for encryption, we must as- 
sess that it is not possible to get an accurate enough estimation of control 
parameters and/or initial conditions. Therefore, a rigorous and concrete the- 
oretical framework is required to quantify the precision of the procedures for 
the estimation of the control parameter and the initial condition of unimodal 
maps from their symbolic sequences. This paper presents this concretion and 
also shows that some of the theorems in [14] are not totally accurate. In this 
sense, those theorems are not only criticized but also rewritten. 



This paper is organized as follows. First of all, Sec. 2 introduces the class of 
maps under study and the main aspects of their symbolic dynamics. Section 3 
remarks the existence of an inner order for the symbolic sequences of a certain 
class of unimodal maps and a relationship between that order and the order of 
the initial conditions employed in their generation. In Sec. 4 the order of the 
symbolic sequences is rewritten in terms of Gray codes and the concept of Gray 
Ordering Number is introduced. After that, Sec. 5 introduces a subclass of the 
class of considered unimodal maps. This subclass of unimodal maps is defined 
in a parametric way, i.e., their dynamics depend on a control parameter. This 
dependency is analyzed by means of the GON. This study will lead to the 
revision and proof of all theorems in [14]. Finally, Sec. 6 summarizes the main 
results of the present work. 
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2 Scenario 



The work described in this paper is focused on a special class of functions. 
This class is denoted by J 7 . A function / belonging to the class T is defined 
in the interval I = [a, b] for a < b and satisfies: 

(1) / is a continuous function in /. 

(2) f(a) = f(b) = a. 

(3) f(x) reaches its maximum value / max < b in the sub- interval [a m , b m ] C I 
so that a m <b m . 

(4) /(/max) < x c and / (/max) > cl, where x c is the middle point of the interval 

In b 1 ip t = am +^ m 
\yj mi f m j , i.e., jj c g 

(5) /(x c ) > x c 

(6) /(x) is an strictly increasing function in [a, a m ] and an strictly decreasing 
function in [b m , b]. 

Xi+i 




Fig. 1. Tent map. 

Although the work in this paper is focused on the class of functions J 7 , it is 
possible to extend it to other class of functions considering the topological 
conjugacy of maps [10, p. 72]. This other class of functions is named T* and 
any / included in J 7 * has the same properties as those in J 7 with the exception 
of properties (3) and (6), since if / is in J 7 *, then it possesses a minimum value 
in [a m , b m ] and is strictly decreasing in [a, a m ] and strictly increasing in [b m , b]. 
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Hereafter, the function f(x) is considered as a way to generate a sequence of 
numbers {x{\ from a certain initial value Xq. Each number Xi determines the 
next element of the sequence trough Xi + \ = f(xi). After a transient number of 
iterations, all the X{ values are inside the interval [x min ,x max ], where x max = 



The tent map is included in the class T and is represented in Fig. 1 . In this case 
a-m = b m = x c and / max = f(x c ) = b. A certain value x i+ i ^ x c can be derived 
from two different values of Xi, as Fig.l informs. In other words, it is satisfied 
that Xi + \ = f(xf) = /(xf), where xf ^ xf , x\ < x c and xf > x c . This is a 
common characteristic of all the functions of the class J 7 . It means that the 
initial condition used in the generation of {x{\ using f(x) can be recovered 
from the last number of the sequence only if the relative position of every Xi 
with respect to x c is known. Therefore, the recovering of the initial condition 
demands recording those relative positions. This is achieved by transforming 
{x^ into a symbolic sequence or pattern according to the next criterion: 

x l = L if xi G [a,x c ), (1) 
Xi = C if Xi = x c , (2) 
x, t = R if Xi G (x c ,b). (3) 

If / is in J 7 * instead of being in J 7 , then the symbolic sequences are generated 
in the same but changing all the L's into i?'s and viceversa. 

Consequently, {xi} is associated to the symbolic sequence P = p pi . . . where 
Pi G {L,R}. Using P and the last element of {xj} one can recover the initial 
condition x . 



3 Relationship between the symbolic sequences and the initial con- 
dition used in their generation 

Let us call Pf(xo) to the symbolic sequence of length n generated from xq 
using the function f(x), which is included in the class T . The value of the 
i— th symbol of the symbolic sequence Pf(xo) is determined by f^(xo), i.e., 
the i— th iteration of f(x) from Xq for i G [0, n — 1]. If pi is the i— th symbol 
of the symbolic sequence, Pi is equal to L if and only if f^(x ) < x c . In the 
same way, pi is equal to R if and only if f^\xo) > x c . As a consequence, the 
definition interval / is divided into 2 t+1 symbolic sub-intervals. Indeed, if x^' 
is the j-th solution of the equation 



f(x c ) and x 





= x, 



(4) 
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the set ja^'-^j f° r < j < 2* divide the definition interval into 2 %+1 sub- 
intervals, where x^, 0,0 -* = x c . All the values included in one of these intervals 
generate the same symbolic sequence of length i + 1. In Fig. 2 the symbolic 
intervals of the tent map for zero, one and two iterations are depicted. The 
main result of the previous proposition is that, for a certain number of itera- 
tions, the different sub-intervals are so that two neighboring sub-intervals lead 
to the same symbolic sequence except for one symbol. On the other hand, for 
i G {0, 1, 2, . . .} and j G [0, 2 % — 1], the set of points x™' determine periodic 
symbolic sequences of period i + 1 when they are considered as initial condi- 
tions. If the symbol C is assigned to x c and only one period is regarded, the 
symbolic sequences generated from jrr^H end with a C. In this sense, if the 
iteration process associated to the generation of a symbolic sequence stops 
just when a C is obtained, only the symbolic sequences derived from the set 
of initial conditions solution of Eq. (4) have finite length. 

All the previous observations can be formally expressed by the following defi- 
nition: 

Definition 1 For a certain function f(x) the symbolic sequence or kneading 
sequence generated from the initial condition xq is Pj(xq). If exists i e No 
such that f^(xo) = x c , then Pf(xo) is finite length. Otherwise, Pf(x$) is a 
kneading sequence of infinite length. As a consequence, any kneading sequence 
of finite length always ends with a C . 

If S is the set of all sequences derived from the iteration of the functions 
included in J 7 , then it is possible to derive a complete ordered set (S, <s) 
where the referred order is defined according to [8, p. 309] as follows: 

Lemma 1 Assuming L <$ C <$ R, S = {s^ and T = {ti} G S, and j is the 
first index so that Sj ^ tj, it is said that S <$ T if one of the next conditions 
is satisfied: 

(1) j = and So <s to- 

(2) j > 0, sqSi . . . Sj-i = toti . . .tj-i contains an even number of R's and 
Sj < s tj. 

(3) j > 0, s si . . . Sj_i = t ^i • • • tj-i contains an odd number of R's and 
Sj > s tj. 

The inner order of S is directly linked to the order on K of the real numbers in 
I used to generate the symbolic sequences from any / in J 7 . This is informed 
and proved in [8, Lemma 4.1] and in [13, Theorem 2]. For the sake of clarity, 
the relationship between the order of the kneading sequences and the order of 
the initial conditions is rewritten as a theorem: 

Theorem 1 For f(x) belonging to the class of functions J 7 and x,y included 
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Fig. 2. Symbolic intervals for different iterations of the tent map. 
in the interval of definition of f(x) so that x < y, it is verified that Pf(x) <s 

Pf(y)- 



4 Gray codes and symbolic sequences 



In the previous section it was remarked that f (jl \x) can be divided into 2 n+1 
intervals such that all the values included in one of those intervals lead to the 
same symbolic sequence of length n + 1. In this sense, those intervals were 
referred as symbolic intervals, since a certain interval can be named through 
the symbolic sequence generated from any value inside it. It was also observed 
that two contiguous symbolic sequences differed in just one symbol. Finally, 
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if the first symbol of the symbolic sequences is discarded, the 2 n+1 symbolic 
sub- intervals generated by the n— th iteration of the map f(x) are symmetric 
with respect to x = x c . In communication theory it is very well known a family 
of codes distinguished by the fact that two successive codes only differ in one 
bit. This family of codes is the Gray codes family, which also presents the 
above cited mirroring property. Table 1 shows the Gray codes of length 4. As 
a result, it is immediate the translation of the symbolic sequences of the class 
of functions J 7 into binary sequences just changing the symbol L into and 
the symbols R and C into 1 [4]. In this sense, the Gray code associated to a 
certain pattern Pf(x) is given by the next definition. 

Definition 2 The Gray code corresponding to Pf(x) = PoPi ■ ■ ■ Pj-i ■ ■ ■ is 
defined as G(Pf(x)) = g^gi ■ ■ ■ gj-i ■ ■ ■ where 



9i 



1 if pi = R 
ifpi = L, 



for i G No- If Pj = C for any j in N , then the Gray code associated to Pf(x) 
is g gi ■■■g j . 

As Table 1 informs, it is possible to translate a Gray code into a binary code. 
The equivalent binary code of a given Gray code can be easily obtained using 
the next definition: 

Definition 3 If the Gray code of a certain symbolic sequence Pf(x) = PqPi ■ ■ -Pj-i 
is given by G(Pf(x)) = g§g\ ■ ■ ■ gj-i • ■ ■ , then the binary code related to Pf(x) 
is U(Pf(x)) = uoU\ ■ ■ ■ Uj-i ■ ■ ■ where 

u i+ i = Ui@ g i+l , 

for i G No and Uo = go. If P/(x) is of length j, i.e., if Pj-i = C, then the 
binary coded related to Pf(x) is U(Pf(x)) = uqU\ ■ ■ -Mj-i where 



u i+ i 



Ui®g i+ i, forQ < i <j - 1, 
1, f ri=j-l. 



Since a binary code can be interpreted as a decimal number just changing the 
base, it is possible to associate a number to a symbolic sequence. However, 
the canonical base changing makes the first symbol modify its weight as the 
length of the symbolic sequence increases. In order to avoid the changing of the 
symbol weights as the length of the symbolic sequences increases, the Gray 
code associated to a symbolic sequence is interpreted as a decimal number 
with integer part equal to zero. The next definition introduces how to carry 
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out the transformation of a symbolic sequence into a real number between 
and 1. 



Rank 


Binary code 
. 


Gray code 





0000 


0000 


1 


0001 


0001 


2 


0010 


0011 


3 


0011 


0010 


4 


0100 


0110 


5 


0101 


0111 


6 


0110 


0101 


7 


0111 


0100 


8 


1000 


1100 


9 


1001 


1101 


10 


1010 


1111 


11 


1011 


1110 


12 


1100 


1010 


13 


1101 


1011 


14 


1110 


1001 


15 


1111 


1000 



Table 1 

Correspondence between Gray codes and binary codes for four bits. 

Definition 4 Let G(P) = g$gi- • • g n -i be a set of bits representing a Gray 
code of length n. Let U(P) = UqUi ■ ■ ■ u n _\ be the binary code corresponding to 
G(P). The Gray Ordering Number or GON of P is defined as the real 
number given by 

GON(P) = 2- 1 ■ u + 2- 2 ■ Ul + ■ ■ ■ + 2~ n ■ 

The definition of the GON also implies the definition of an order <gon upon 
the set of symbolic sequences S. In other words, according to the definition 
of the GON, it is possible to build the complete ordered set (<S, <gon)- This 
ordered set is equivalent to (S, <s), i-e., the order defined using the GON is 
equivalent to the order < 5 . 

Proposition 1 The orders <$ and <gon (ire equivalent on S. 
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PROOF. Let P = poPi ■ ■ .pj-ipj ... be a certain symbolic sequence which 
can be of finite or infinite length. If U (P) = UqUi . . . Uj^iUj ... is the binary 
code linked to the kneading sequence P, Uj is equal to 1 if one of the next 
situations occurs: 

(1) pj = R and poPi ■ ■ -Pj-i contains an even number of i?'s. 

(2) pj = L and poPi ■ ■ -Pj-i contains an odd number of i?'s. 

Let Q = q^qi . . . qu-iQk • • • be another kneading sequence of finite or infinite 
length. Let U(Q) = t^ti . . . t^-itk ... be its associated binary code. According 
to Theorem 1, if the first different symbol between P and Q is the i— th one, 
then P <s Q if and only if one of the next cases happens: 

(1) pi = R, qi = L and p pi . . -Pi-i contains an odd number of i?'s. As a 
consequence, it is verified that = and ti = 1, which implies that 
GON(P) < GON(Q), i.e., P < GON Q- 

(2) pi = R, Q of length i and p pi . . . contains an odd number of i?'s. 
Since Q is finite-length, its final symbol is C. Therefore, t« = 1 and Ui = 
implying that GON(P) < GON(Q), i.e., P < GON Q- 

(3) pi = L, qi = R and poPi • • • Pi-i contains an even number of i?'s. For this 
configuration, Ui = and = 1, which informs GON(P) < GON(Q) 
and subsequently P <gon Q- 

(4) P of length i, qi_ x = R and p pi . . -Pi-2 contains an even number of i?'s. 
Since P has % symbols, it means = 1. On the other hand, = 1 
and three possible situations are posible 

(a) Q is of length j for j > i. Then = 1 implies GON(P) < GON(Q. 

(b) Q is infinite-length and g« = L, implying £j = 1 and GON(P) < 
GON(Q). 

(c) Q is infinite-length and g« = R. In this case there exists j > % such 
that qj = R. Otherwise, the condition P <$ Q implies that P is of 
length 1 and Q = RLLLL .... In each of these situations it is satisfied 
GON(P) < GO N {Q). 

On the other hand, let us assume P <gon Q and i the first index such that 
Ui ^ t { . 

(1) Ui — 0, ti — 1 and popi . . -Pi-i contains an odd number of i?'s. Since 
GON(P) < GON(Q), then pi = R and = L, which further implies 
that P < S Q. 

(2) Ui — 0, ti — 1 and p^pi . . - Pi-i contains an even number of i?'s. In this 
situation the assumption GON(P) < GON(Q) forces Pi = L and g« =r, 
which informs that P <$ Q. 

(3) P is of length i, t^i = 1. This implies that p^pi . . - Pi-2 contains an even 
number of i?'s, q^i = R and thus P <s Q- 

(4) Ui-i = and Q of length i and p Pi • • - Pi-i contains an odd number of 
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i?'s. Therefore, = R, g.j_i = C and P <s Q- 
As a result, P <s Q if and only if P <gon Q and the proof is complete. □ 

The previous proposition and Theorem 1 lead to the next theorem, which 
represents the extension and proof of Theorem 1 in [14]. 

Theorem 2 For f G 7 and x,y G I, it is satisfied that GON(Pf(x)) < 
GON(Pf(y)) if and only if x < y. In other words, the GON of the symbolic 
sequences in S is an increasing function with respect to the initial condition. 



5 Gray codes and parametric unimodal maps 

A special case of interest is the study of unimodal maps defined in a parametric 
way. In this sense, this section is focused on the analysis of the class of functions 
fx(x) G 7 for all A in [0,1]. Let F(x) G 7 and F(x c ) = F max < b. The 
parametric function fx can be expressed as follows: 

fx(x) = XF(x), (5) 

which implies fx(x c ) = X ■ F max , which is the maximum value of fx(x). A first 
consequence of this is Theorem 3 in [14], which is a corollary of Theorem 2. 

Corollary 2.1 For f x (x) = \F(x) with F(x) G J 7 and A G [0, 1], it is satisfied 
that GON{P h {f x {x))) < GON{P h {fx{x c ))),\/x G [a,b\. 

Moreover, the maximum value of fx(x), i-e., AF max depends on A in such 
a way that an increment of the control parameter forces an increment of 
the maximum value. As a consequence, the GON of the kneading sequences 
derived from x = fx{x c ) is an increasing function with respect to the control 
parameter [J I, Theorem 4]. 

Corollary 2.2 For f x (x) = XF(x) with F(x) G 7 and X U X 2 G [0,1] with 
Aj < X 2 , it is satisfied that GON(P hi (f Xl (x c ))) < GON{P fx2 {fx 2 {x c ))). 

On the other hand, after a certain number of transient iterations, all the values 
obtained from any initial condition through the iteration of any function in 7 
are inside the interval [x min , x max ]. Therefore, once all the values derived from 
the iteration of the considered function are inside [x min , x max ], it is verified 

(2) 

that GON(P h (x)) > GON(P h (fX>(x c )). This was wrongly interpreted in 
[14, Theorem 5], since this theorem is only satisfied if f^\x) > x m i n for any 
x G [a,b]. Nevertheless, the previous comments point out that this inequality 
is verified only for x G [/^ 1 (/^ 1 (^min)), b], i-e., Theorem 5 in [14] is not fulfilled 
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for x G [ a , f\ {fx ( x mm)))- Consequently, it is necessary to modify Theorem 
5 in [] : ] according to the preceding considerations. In this sense, the next 
corollary rewrites Theorem 5 in [1 i] in a more accurate way and, at the same 
time, extends its application domain to all the functions in F . 

Corollary 2.3 Let F(x) be a function in F that leads to f\(x) = XF(x) for 
x G [a, b] and A G [0, 1]. Let Xi be defined as Xi = x for i = and X{ = f\(xi-i) 
for i > 0, i G N. There exists n\ G N such that Xi is in [x m i n , x max ] for i > n\ 
and it is satisfied that GON(P h (xi)) > GON(P fx (fj® (x c )),Wx G [a, b] for 
i > n\. 

(2) 

Finally, the value x min is given by frixc) = h(fx(x c )) = A • F(AF max ). If 
x m i n is a monotonic function of A, then it is possible to extract a new corollary 
from Theorem 2. In [1 1, Theorem 6] it is assumed without proof that f\\x c ) 
is a monotonic decreasing function with respect to A. This assumption implies 
that 



F(AF„») + A.F„ 9F[X) 



OA OX 



< 0. (6) 



ax 

This condition is not satisfied for all the possible values A and for all the 
functions in T . Let us consider the logistic map. In [14] the dependency of 
x min on A is studied using the logistic map. Indeed, the logistic map is a 
function included in J 7 , which is defined as 

f x (x) = X-Ax(l-x), (7) 

for A G [0, 1] and x G [0, 1]. It is easy to verify that for the logistic map the 
condition given by Eq. (6) is fulfilled if and only if A > 8/12. Therefore, Theo- 
rem 6 in [14] must be rewritten in such a way that the discussed inaccuracy is 
overcome and, simultaneously, the application domain of its variant affects not 
only the logistic map but all the functions in J 7 . Again, this aim is completed 
through a series of additional assumptions on the scope defined in Theorem 
2. 

Corollary 2.4 Let us suppose that f\(x) = XF(x) with F(x) G F , X G [0, 1] 

(2) 

and x G [a, b\. For Ai, A2 G [0, 1] with Ai < A2 and satisfying df{ (x c )/dX < 
forX = {X X ,X 2 }, itis verified that GON(P hi (f^(x c ))) > GON(P fx2 (f { *>(x c ))). 



6 Conclusions 



In this paper we have mathematically proven that it is possible to read the 
classical theory of applied symbolic dynamics for unimodal maps from the 
point of view derived from the concept of Gray Ordering Number. Indeed, 
the main results of the present work were previously presented in other works 
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as theorems. Nevertheless, these theorems were not formally demonstrated. 
We have provided not only the mathematical proof of these theorems but 
also solved some imprecisions, which is essential to use the concept of Gray 
Ordering Number in a correct and efficient way. The main result of all this work 
is the possibility of improving and expanding previous contributions based on 
the concept of Gray Ordering Number. Specially relevant is the case of the 
estimation of the values of the initial condition and the control parameter of 
unimodal maps. The theoretical framework presented in this paper allows to 
establish the limitations of the methods previously proposed for the estimation 
of those values. Furthermore, this paper is the theoretical conclusion of all the 
work that we have carried out on unimodal maps both in the field of the 
applied theory of symbolic dynamics [4,3,5], and in the context of chaos-based 
cryptography [2,7,12,6]. 
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